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Abstract 

Mean field dynamo theory is a leading candidate to explain the large scale 
magnetic flux in galaxies and stars. However, controversy arises over the ex- 
tent of premature quenching by the backreaction of the growing field. We 
distinguish between rapid dynamo action, which is required by astrophysical 
systems, and resistively limited dynamo action. We show how the flow of 
magnetic helicity is important for rapid dynamo action. Existing numeri- 
cal and analytic work suggesting that mean field dynamos are prematurely 
quenched and resistively limited include approximations or boundary condi- 
tions which suppress the magnetic helicity flow from the outset. Thus they 
do not unambiguously reveal whether flux generating astrophysical dynamos 
are dynamically suppressed when the required helicity flow is allowed. An 
outflow of helicity also implies an outflow of magnetic energy and so active 
coronae and/or winds are a prediction of mean field dynamos. Open bound- 
aries alone may not be sufficient and the additional physics of buoyancy and 
winds may be required. Some possible simulation approaches, even with 
periodic boxes, to test the principles are discussed. Some imitations of the 
"Zeldovich relation" are also addressed. 
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I. INTRODUCTION 



A leading candidate to explain the origin of large scale magnetic flux 
growth in stars and galaxies has been the mean field dynamo (MFD) the- 
ory [1,2,3,4,5,6]. The theory appeals to a combination of helical turbulence 
(leading to the a effect), differential rotation (the Q effect), and turbulent dif- 
fusion to exponentiate an initial seed mean magnetic field. Ref [7] developed 
a formalism for describing the concept [8] that helical turbulence can twist 
toroidal (0) fields into the poloidal (r, z) direction, where they can be acted 
upon by differential rotation to regenerate a powerful large scale toroidal 
magnetic field. Their formalism involved breaking the total magnetic field 
into a mean component B and a fluctuating component b, and similarly for 
the velocity field V. The mean can be a spatial mean or an ensemble av- 
erage. For comparison to observations of a single astrophysical system, the 
ensemble average is approximately equal to the spatial average when there 
is a scale separation between the mean scale and the fluctuating scale. In 
reality, the scale separation is often weaker than the dynamo theorist desires. 
Nevertheless, we proceed to consider spatial averages. 

Ref [7] showed that B satisfies the induction equation 

'YD 

^ = - C VxE, (1) 

where 

E = - (V x B") /c - (v x b)/c + AV x B, (2) 
where the first term describes the effect of differential rotation ("fi- effect"), 

(v x b)i = OijBj - PijkdjBk (3) 

is the "turbulent emf," and A = r/c 2 /4rr is the magnetic diffusivity defined 
with the resistivity 77. Here contains Parker's twisting ("a effect") and 
Pijki^ A) contains the turbulent diffusivity. Ref [7] calculated E to first 
order in B for isotropic o;^ and flijk and hence the dynamo coefficients a and 
(3 to zeroth order in B from the statistics of the turbulence. They ignore 
the Navier-Stokes equation. The back-reaction on the dynamo coefficients to 
first order in B and V was calculated in Ref [9], and Ref [10] calculates a 
to all orders in B when mean field gradients are small) See also Refs [11, 12] 
and the issues raised in [13]. Refs [11,12] obtain a catastrophically quenched 
a but by a different argument than we use later. 
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When (0) is substituted into ([I]), we have the mean-field dynamo equation: 

'YD 

_ = Vx(VxB)+Vx («B) - V x (p + A) V x B. (4) 

In the approximation that V, a and (3 are independent of B, (|4|) is a linear 
equation for B which can be solved as an eigenvalue problem for the growing 
modes in the Sun and other bodies. Actually, a rapid growth of the fluctu- 
ating field necessarily accompanies the mean-field dynamo. Its impact upon 
the growth of the mean field, and the impact of the mean field itself on its 
own growth are controversial. 

The controversy results because Lorentz forces from the growing magnetic 
field react back on and complicate the turbulent motions driving the field 
growth [14,15,16,17,18,19]. It is tricky to disentangle the back reaction of the 
mean field from that of the fluctuating field. Analytic studies and numerical 
simulations seem to disagree as to the extent to which the dynamo coefficients 
are suppressed by the back reaction of the mean field. But, as we will address, 
the disagreements may be because different problems are being solved and 
the results inappropriately compared. 

It is important to distinguish between rapid MFD action, resistively lim- 
ited MFD action, and no MFD action. Rapid dynamo action describes dy- 
namo action which proceeds at rates not greatly suppressed from the kine- 
matic values. Resistively limited dynamo action is MFD growth that pro- 
ceeds at rates strongly dependent on the magnetic Reynolds number. No 
MFD action implies no mean growth at all. For galaxies, rapid MFD action 
is necessary if the observed large scale fields are to be produced and sustained 
by the MFD. There the a effect and the Q effect operate in the same volume. 
For the Sun, rapid dynamo action is also required but the a effect operates 
in the convection zone, whilst the Q shear effect operates in the overshoot 
layer beneath the convection zone [20]. 

In section II we discuss the role of magnetic helicity in dynamo theory. 
We show that open boundaries, allowing magnetic helicity to escape, play 
an important role for rapid astrophysical dynamo action. This leads to an 
ambiguity in interpreting of quenching simulations. In section III we predict 
that astrophysical rotators with dynamos harbor steady coronae. In section 
IV we suggest that open boundary conditions may not be enough for a dy- 
namo and that the magnetic helicity flow may require buoyancy or winds. 
We also address pitfalls of the Zeldovich relation. We conclude in section V. 
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II. ROLE OF MAGNETIC HELICITY CONSERVATION 



Although the MFD theory predates detailed studies of MHD turbulence, 
the MFD may be looked upon in hindsight as a framework for studying the 
inverse cascade of magnetic helicity. Whether this inverse cascade is primarily 
local (proceeding by interactions of eddies/waves of nearby wavenumbers) 
or non-local (proceeding with a direct conversion of power from large to 
small wavenumbers) is important to understand. The simple MFD is most 
consistent with the latter. 

From the numerical solution of approximate equations describing the 
spectra of energy and helicity in MHD turbulence, Ref [21] showed that the 
a effect conserves magnetic helicity (= J (A - B)d 3 x) by pumping a positive 
(negative) amount to scales > L (the outer scale of the turbulence) while 
pumping a negative (positive) amount to scales <C L. Magnetic energy at 
the large scale was identified with the B of [7]. Thus, dynamo action leading 
to an ever larger B, hence the creation of ever more large scale helicity, can 
proceed as long as small scale helicity of opposite sign can be removed or 
dissipated. Recent simulations [22] confirm this inverse cascade and the role 
of H M conservation. But the rate of small scale H M removal determines the 
rate of MFD action. Presently, simulations have invoked boundary condi- 
tions for which the growth of large scale field is resisitvely limited. We now 
discuss what this means. 

A. Constraint equation for the turbulent EMF 

To explore the role of boundary conditions in constraining the value of 
the a dynamo parameter, we take Ohm's law 

E=-c 4 VxB + ? )J (5) 

and average the dot product with B to find [23] 

(E-B) = E-B + (e-b) = -c _1 (v x b) • B + 77J • B + (e • b) (6) 

where J is the current density. 

A second expression for (E • B) also follows from Ohm's law without first 
splitting into mean and fluctuating components, that is 

(E • B) = 77 ( J • B) = r]3 • B + 77 (j • b) = 77T ■ B + c -1 A (b • V x b) . (7) 
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Using (0) and (||), we have 

-c 4 (vxb)-B = c^X (b • V x b) - (e ■ b), (8) 

which can be used to constrain (v x b) in the mean field theory. 

B. Necessity of magnetic helicity escape and open boundaries 

Now consider E in terms of the vector and scalar potentials A and $: 

E = -V$ - (l/c)<9 t A. (9) 

Dotting with B = V x A we have 

E-B = -W-B-(l/c)B-d t A. (10) 

After straightforward algebraic manipulation, application of Maxwell's equa- 
tions and V • B = 0, this equation implies 

E B = -(1/2)V • $B + (1/2)V • (A x E) 
-(l/2c)3t(A • B) = (-l/2c)«V^(B) = V 3 B, (11) 

where 

H^B) = (H , H) = [A ■ B, c$B - cA x E] (12) 

is the magnetic helicity density 4- vector [24], and the contraction has been 
done with the 4x4 matrix r]^ u where r]^ = for /i ^ u, r]^ — 1 for // = v = 
and 77 M „ = — lfor/i = z/>0. Taking the average of ([□]) gives 

9 M F M (B) = -2c(E ■ B) = -2cE • B - 2c(e • b) = -2cr](J ■ B). (13) 

If, instead of starting with the total E as in ([|), we start with e and then 
dot with b and average, the analagous derivation replaces (|TJ) by 

d v W{h) = -2c(e- b), (14) 

where H^ih) indicates the average of H^(h). The latter is defined like (|12"D 
but with the corresponding fluctuating quantities replacing the total quanti- 
ties. Similarly, starting with E and dotting with B, gives 

df J ,H tl (B) = — 2cE • B, (15) 
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where if ^(B) is defined as in fli"2| ) but with the corresponding mean quantities 
replacing the total quantities. Now consider two cases. 

1. Case 1 



In this case we consider statistically stationary turbulence in which the 
scale of averaging is equal to the universal scale, or equivalently for present 
purposes, that the overbar indicates averaging over periodic boundaries. In 
this case, the spatial divergence terms on the left of (|T3"D become surface 
integrals and vanish. Rewriting this ([14]), and discarding the divergenece 
terms gives 

<e-b) = -l«9 t (a-b> (16) 

Thus, eqn (|j) is resistively limited unless significant non-stationarity is al- 
lowed. 

Let us apply this to the specific case of statisically stationary turbulence 
in which a uniform mean field is imposed over a periodic box, and for which 
the averaging scale is the scale of the box. Then the mean field cannot 
change with time, and has no gradients. This is the case of Ref [18]. Then 
the right side of (|TB]) vanishes because not only are zeroth order turbulent 
correlations (those computed to zeroth order in the mean field) stationary, 
but all higher order corrections must also be stationary. The mean field does 
not change with time so no mean quantity varys on macroscopic time scales. 
Incidentally, when expanded in powers of the mean field, the lowest order 
contribution to (|T6|) enters at second order [23]. 

For the uniform field in a periodic box, (§) then implies that 



- c _1 (v x b) ■ B = aW/c = c _1 A (b • V x b) 



(17) 



where a = for a uniform field in the z direction. Rearranging we have 

A (b ■ V x b) 



w 



(18) 



Now if we assume that b and v follow e.g. a k 5//3 Kolmogorov energy spec- 



trum, we then have 



ft < 



k blr 




ft 


<B> 2 







(19) 
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where n = 3/4 if the current helicity is dominated by large wavenumber and 
n = 1 if it is dominated by small wavenumber. The v , b are the speed and 
fluctuating field magnetic energy of the dominant energy containg eddies, ko 
is the wavenumber for that scale, tq is the associated eddy turnover time 
and the magnetic Reynolds number Rm = vL/X, where L is the scale of 
the energy dominating eddies. The latter similarity in follows for near 
equipartition between magnetic and kinetic energies on the outer scale, and 
we have written the kinematic a as o.q. (Note that the "Pouquet correction" 
[21] to a which is oc (b • V x b) enters to the same order in the mean 
field as the standard (v • V x v) term [10,32], so we consider either to be 
representative of the lowest order contribution). The value of a in this case 
is "resistively limited," however the Reynolds number factor on the bottom 
does not represent a dynamical backreaction. It is an a priori implication of 
the imposed boundary conditions. Thus simulations which invoke periodic 
boundary conditions and find this level of suppression [18], may not be seeing 
the effect of backreaction dynamics, but of the boundary conditions. 

2. Case 2 

In this case, consider the system (e.g. Galaxy or Sun) volume V << 
universal volume. Integrating ( JTT| ) over all of space, U, then gives 

f E • B d 3 x = -(1/2) J V • $B d 3 x + (1/2) J V • (A x E) d 3 x 

-(l/2c)d t [ A ■ B d 3 x = -(l/2c) W(B) = / riJ ■ Bd 3 x, (20) 
Ju Ju 

where the divergence integrals vanish when converted to surface integrals at 
infinity. We have defined the global magnetic helicity 



ft(B) = / A • B d 3 x, (21) 



iu 

where U allows for scales much larger than the mean field scales. It is straight- 
forward to show that a parallel argument for the mean and fluctuating fields 
respectively leads to 



d t K(B) = 9tj A- B d 3 x = -2c^E-B d 3 x (22) 



and 



dt'Kih) = d t f (a • b) d 3 x = -2c / (e ■ b) d 3 x = -2c / e • b d 3 x = cWb), 
Ju Ju ' Ju 

(23) 
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where the penultimate equality in (^) follows from redundancy of averages. 

We now split ( j2"2"|) and (|23|) into contributions from inside and outside the 
rotator. One must exercise caution in doing so because Oi, is gauge invariant 
and physically meaningful only if the volume U over which J-C is integrated 
is bounded by a magnetic surface (i.e. normal component of B vanishes at 
the surface) , whereas the surface separating the outside from the inside of 
the rotator is not magnetic in general. 

Ref [25] shows how to construct a revised gauge invariant quantity called 
the relative magnetic helicity. This can be written as 

ftii^Bi) = W(Bi, P ) - ft(Pi, P ) (24) 

where the two arguments represent inside and outside the body respectively, 
and P indicates a potential field. The relative helicity of the inside region 
is thus the difference between the actual helicity and the helicity associated 
with a potential field inside that boundary. The use of P; is not arbitrary 
in (|24j ), and is in fact the field configuration of lowest energy. While (j2jj) is 
insensitive to the choice of external field [25], it is most convenient to take 
it to be a potential field as is done in (^) symbolized by P D . The relative 
helicity of the outer region, !Kr i0 , is of the form (|24] ) but with the o's and i's 
reversed. The "K R is invariant even if the boundary is not a magnetic surface. 

The total global helicity, in a magnetically bounded volume divided into 
the sum of internal and external regions, U = Ui + U e , satisfies [25] 

K(B) = K Rj0 (B)+K R>i (B), (25) 

when the boundary surfaces are planar or spherical. This latter statement 
on the boundaries means the vanishing of an additional term associated with 
potential fields that would appear in (p5|) . Similar equations apply for B and 
b, so (|22| ) and ( p3|) can be written 

c^B) = d t H R , (B) + d t K R>i (B), (26) 

and 

d t %{b) = d t H R , {b) + d t M R>i {b) (27) 
respectively. According to equation (62) of Ref [25], 

d t 3t Rji (B) — -2c / B-Bd 3 x + 2cf (A p x E) ■ dS, (28) 
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where A p is the vector potential corresponding to a potential field P in 
U e , and DUi indicates integration on the boundary surface of the rotator. 
Similarly, we have 

dt'KRiCB) = -2c [ E ■ Bd 3 x + 2c [ (A„ x E) • dS (29) 
JUi Jdu % 

and 

d t K Ri (b) = -2c I e-hd 3 x + 2c[ a^xe-rfS. (30) 

JU, JDU t 

Note again that the above internal relative helicity time derivatives are both 
gauge invariant and independent of the field assumed in the external region. 
If we were considering the relative helicity of the external region, that would 
be independent of the actual field in the internal region. 

Now if we take the averaging scale to be less than or equal to the Ui scale, 
we can then replace (|30|) by 

(e ■ b) = -l^JC^b) + (V • (a p x e)>, (31) 

where the brackets indicate integrating over Ui or smaller. We now see that 
even if the first term on the right of ([JI]) vanishes, (e • b) contributes a surface 
term to (H) that need not vanish. The turbulent EMF is not resistively limited 
as in the case of the previous section. The surface term can strongly dominate 
the resistive contribution. Thus in a steady state for Rm » 1, an outflow 
of magnetic helicity is likely essential to keeping the dynamo operating. For 
the Sun, the steady state would refer to time scales longer than the longest 
eddy turnover time, but shorter than the predicted 11 year Solar cycle. Note 
that the right of ([13]) is small for for large magnetic Reynolds numbers, so 
the flux of helicity has contributions from the small and large scale field. If 
the surface term vanishes, helicity could instead be "injected" through the 
time derivative term of (|13D. This has known application to Tokomaks [26]. 

In short, steady dynamo action unrestricted by resistivity is possible in 
case 2 but not in case 1. This is consistent with current simulations. Case 

1 is represented by Refs [21] and [22], which find that when the boundary 
conditions are periodic, dynamo action proceeds at resistively limited rates. 
In [27] action seems to proceed more rapidly, and this corresponds to a case 

2 simulation, though the resolution and Rm were small. 

III. ENERGY FLOW TO CORONAE 
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Here we explore only the deposition of relative helicity to the exterior 
and the associated total magnetic energy without addressing how the energy 
converts to particles or flows. We assume that the rotator is in a steady state 
over the time scale of interest, so the left sides of (p9|) and (|30| ) vanish. Note 
that in a system like the Sun where the mean field flips sign every ~ 11 years, 
the steady state is relevant for time scales less than this period, but greater 
than the eddy turnover time (~ 5 x 10 4 sec). Beyond the ~ 11 year times 
scales, the mean large and small scale relative helicity contributions need not 
separately be steady and the left hand sides need not vanish. 

The helicity supply rate, represented by the volume integrals (second 
terms of ( p9|) and (|30|) ), are then equal to the integrated flux of relative 
magnetic helicity through the surface of the rotator. Moreover, from (|T3"D, 
we see that the integrated flux of the large scale relative helicity, = 3^(13), 
and the integrated flux of small scale relative helicity, = ^^(b), are equal 
and opposite. We thus have 

9 r R)i (B) = -2c f E • Bd 3 x = 2c f e • bd 3 x = -^(b). (32) 

JUi J Ui 

To evaluate this, we use ([[]) and (0) to find 

E= -c -1 (aB-/3V x B), (33) 

throughout Ui. Thus 

5F^(B) = -9^(b) = 2 / («B 2 - fffi ■ V x B)d 3 x. (34) 

This shows the relation between the equal and opposite large and small scale 
relative helicity deposition rates and the dynamo coefficients. 

Now the realizability of a helical magnetic field requires its turbulent 
energy spectrum, E^ 1 , to satisfy [28] 

£f (b) > ±k\X k Q>)\ , (35) 

where %k is the magnetic helicity at wavenumber k. The same argument 
also applies to the mean field energy spectrum, so that 

£f (B) > ^-k\0i k (B)\. (36) 
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If we assume that the time and spatial dependences are separable in both 
E M and IK, then a minimum power delivered to the corona can be derived. 
For the contribution from the small scale field, we have 



E M (h) = JEf(b)dk> ±fk\? ktR ,i(b)\dk> ^j|%(b)|_* 

>%HJ^(b)| = %HJ^(B) 



where the last equality follows from the first equation in (|34]). The last quan- 
tity is exactly the lower limit on E M (B). Thus the sum of the lower limits 
on the total power delivered from large and small scales is 13^(13)1 + 



''''" >n |IF Rii (b)| = 2^|3 r (B)|. Now for a mode to fit in the rotator, k > k min = 
27r//i, where h is a characteristic scale height of the turbulent layer. Using 
(|34|), the total estimated energy delivered to the corona (=the sum of the 
equal small and large scale contributions) is then 

EM> 2 ^\? (b) | =2 *p| ? (b)| = L B 2 -0-VxB , (37) 

87T 87T n ave 

where V is the volume of the turbulent rotator and ave indicates volume 
averaged. We will assume that the two terms on the right of Q57D do not 
cancel, and use the first term of (|37j ) as representative. 

Working in this allowed time range for the Sun (5 x 10 4 sec < t < llyr 
we apply (j57p to each hemisphere of the Sun, and using the first term as an 
order of magnitude estimate, obtain 

(3S) 

where we have taken a ~ 40 cm s _1 , and we have presumed a field of 150G 
at a depth of 10 4 km beneath the solar surface in the convection zone, which 
is in energy equipartition with turbulent kinetic motions [2]. 

As this energy deposition rate is available for reconnection which can 
generate Alfven waves, drive winds, and energize particles, we must compare 
this limit with the total of downward heat conduction loss, radiative loss, and 
solar wind energy flux in coronal holes, which cover ~ 1/2 the area of the 
Sun. According to [29] this amounts to an approximately steady activity of 
2.5 x 10 28 erg s~ l , about 3 times the predicted value of (p8|) . Other supporting 
evidence for deposition of magnetic energy and magnetic + current helicity 
[30] in the Sun is discussed in Ref [31]. 
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AGN and the Galactic ISM represent other likely sites of mean field dy- 
namos. Ref [31] discusses the associated energy deposition to the coronae of 
these systems as well. For the Galaxy, E M > (nR 2 )aB 2 ~ 10 40 ( J R/12kpc) 2 
x (a/10 5 cm/s)(i?/5 x 10~ 6 G) 2 erg/s in each hemisphere. This is consistent 
with coronal energy input rates required by [32] and [33] . For AGN accretion 
disks, the deposition rate seems to be consistent with what is required from 
X-ray observations. Independent of the above, the most successful paradigm 
for X-ray luminosity in AGN is coronal dissipation of magnetic energy [34]. 

IV. DISCUSSION AND OPEN QUESTIONS 

A. Role of boundary diffusion and winds 

We have seen that a steady MFD unlimited by resistivity requires the he- 
licity to flow through the rotator boundary. However, merely open boundary 
conditions may not be enough: some mechanism must sustain the diffusion 
of the mean field at the boundary. Near the surface, buoyancy or winds [35] 
and winds may have to play a role. Note however, that it is the mean field 
which needs to diffuse, not necessarily the total field, and not the matter [5]. 

For the Galaxy, turbulent diffusion of the mean magnetic field across the 
boundary is also required to maintain a quadrupole geometry with a net flux 
inside the disk. Ref [36] argues that surface diffusion of total magnetic field 
for the Galaxy is difficult. Thus the diffusion of the mean and/or total field 
out of a galaxy remains an open question. For the Sun, the solar cycle also 
requires diffusion of mean field through the boundary. The flow of helicity 
would appeal to the same dynamics needed by these constraints. 

The role of boundary diffusion is essential for understanding whether the 
present Galactic field is primordial or produced in situ and how the solar 
dynamo accounts for the surface fields. 

B. Choosing simulation boundary conditions: transient vs steady 
phases 

Ideally, the best simulation to test the dynamo would be a fully global 
simulation with significant scale separation between the system size and the 
turbulent scale, and would properly include the boundary diffusion mecha- 
nisms with open boundary conditions. However it may be possible to invoke 
periodic boundary conditions to simulate rapid dynamo action in the follow- 
ing ways. The first represent transient phases, the last approaches represent 
steady action. 
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A first transient phase during which rapid dynamo action should incur, 
is the time period during which the helicity on the small scale builds up. 
During this period the dynamo can grow rapidly. After this period one is 
subject to the boundary conditions. This has been investigated numerically 
in Ref [37], where the dynamo a is shown not to be significantly quenched 
in this early phase. 

A second time scale over which rapid dynamo action could be seen in a 
periodic box, even after the small scale helicity saturates, would be facilitated 
by significant scale separation between the box size and the scale of the input 
turbulent scale. One could inject kinetic helicity only into a sub-volume of 
the box and observe the inverse cascade of magnetic helicity and magnetic 
energy in this sub-volume before the material diffuses all the way across the 
box. During this transient period, even if the small scale helicity saturates, 
the surface terms can play an important role. The length of this period of 
course depends on the extent of scale separation. 

Finally, there may be a way to use a periodic box to actually get steady 
rapid dynamo action. One can inject helicity inside the sub-volume described 
above, and invoke a one-way diffusive valve out of this region. The helicity of 
one sign should grow inside the forced region, while the opposite sign would 
be shed to the exterior region. One could then surround the box with a 
resistive layer so that this external helicity dissipates rather than feeds back 
into into the system. 

C. Previous analytical and numerical suppression results 

To date, analytic and numerical studies that suggest resistively limited 
dynamo action either (1) invoke periodic boundary conditions (as described 
above), and/or (2) are 2-D, or (3) do not distinguish between zeroth or- 
der isotropic components of the turbulence and the higher order anisotropic 
perturbations for a weak mean field [31]. This existence of an alternative 
explanation for the catastrophic suppression makes the computed suppres- 
sion ambiguous. This does not mean that the physical concepts found in the 
strong suppression results are invalid, but just that they may be valid only 
for the restricted cases considered. 

The limitations can be specifically pinpointed. For example, [38] derives 
the now famous "Zeldovich relation" (b 2 )/B ~ Rm- But this relation is for 
2 — D, ignores boundary terms, and is non-trivial only when the magnetic 
energy is dominated on small scales (unlike 3-D simulation results). 
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The Zeldovich relation arises from first deriving, in the absence of bound- 
ary terms, the evolution equation for the average of the square of the vector 
potential from an initially uniform magnetic field in 2-D. We choose the vec- 
tor potential to be in the z direction, and thus we have A = A z . The resulting 
evolution equation is 

d t (A 2 ) = -2\{{VA Z ) 2 ). (39) 

Note that there is only a single term on the right hand side because of the 
boundary conditions. We now appeal to the positive definite nature of the 
last term to rewrite the equation 

d t (A z 2 ) = -2\(A z 2 )/5\t), (40) 

where 5(t) represents the dominant scale of the turbulent field. In 2-D, when 
surface terms are ignored, the mean field cannot change. However the small 
scale field can grow in response to the turbulent stretching. In analytic 
studies [eg. 16] and numerical simulations [e.g. 37] (which happens to be 
3-D) this type of behavior can be described by a decrease of 6(t), that is 
d5(t)/dt < 0, as the turbulent eddies cascade to smaller and smaller scales. 
Equipartition between magnetic and kinetic energy is first reached on the 
dissipative scale as there the eddy turnover time is the shortest. When the 
field peaks at this scale, the dissipation rate in ( fl"0|) is maximized because 
S(t) reaches its minimum. The end of kinematic regime occurs when the field 
saturates at the resistive scale. In fact the maximum dissipation rate occurs 
if Rm ~ b 2 /B . (Note however, that if R M » b 2 /B , saturation at b 2 = v 2 
would occur before Rm enters the relation.) 

Soon after the maximum dissipation is achieved, d5/dt > 0. That is, the 
scale 5(t) increases as the dynamic regime sets in and equipartition is reached 
on successively larger and larger scales. The dissipation rate then decreases 
in the non-linear dynamic regime. This is the "suppression" of dissipation 
that is directly implied by (0). We refer to this as "dissipation" rather than 
"diffusion" because the only term determining the time evolution of (A 2 ) has 
an explcit A coefficient. 

In fact, the original Zeldovich result was meant to be a limitation on the 
growth of (b 2 ) from a fixed B, the latter of which was restricted not to change 
due to 2-D calculation and the choice of boundary conditions. Nothing was 
said about the turbulent diffusion of field lines or the ability of the mean field 
to grow. Zeldovich recognized that the 2-D case did not necessarily imply 
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much about the 3-D case, the case when boundary terms are included, or the 
case when B is allowed to grow. 

That being said, there are some important related issues to understand in 
future work. The observation that at least the 2-D Lagrangian chaos prop- 
erties of the flow seem to change in the presence of a weak mean field for 
turbulence in a periodic box [39] needs to be understood in relation to the 
imposed boundary conditions and the shape of the magnetic energy spec- 
trum. The connection between this observed restricted 2-D diffusion and the 
Zeldovich relation is actually subtle, and also requires further investigation. 
One point to consider in this context is that when the diffusion time across 
the box is comparable to the Alfven crossing time across the box, numerical 
effects may be playing a strong role in determining the results and must be 
kept in mind. 

At present there is a dearth of simulations on 3-D turbulent diffusion of 
the mean field. However molecular cloud simulations [40] which do not ex- 
plicitly consider an application to dynamo theory, can actually be interpreted 
to imply the absence of catastrophic quenching of mean field diffusion. 

Finally, we note that Zeldovich et al. produced a second approach to this 
relation between the mean field strength and Rm which was thought to be 
applicable to 3-D [4]. In 3-D, a logarithmic dependence on Rm arises, while 
the method also seems to reproduce the 2-D result in the 2-D limit. However, 
this approach may be flawed as it does not distinguish between isotropic and 
anisotropic components of the turbulence. This will be discussed further 
elsewhere. 

D. Coronal Activity 

The estimated energy deposition rates are consistent with the coronal + 
wind power from the Sun, Galaxy and Seyfert Is [31]. The helical properties 
also seem to agree well in the solar case where they can be observed [30]. 
The steady flow of magnetic energy into coronae thus provides an interesting 
connection between mean field dynamos and coronal dissipation paradigms 
in a range of sources. A reasonably steady (over time scales long compared to 
turbulent turnover time scales), active corona with multi-scale helical struc- 
tures, provides a self-consistency check for a dynamo production of magnetic 
field in which there is exponential field growth inside the body on small and 
large scales. If the growth rate were only linear, the corona and wind output 
might be more episodic with fluctutuations not by many orders of magni- 
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tude but by an order of magnitude or less. The latter seems to be consistent 
with observations of the Sun and coronae of Seyfert AGN [41], suggesting 
exponential field growth inside the rotator. 

E. Role of magnetic Prandtl number 

Much of the numerical and analytic work to date with respect to the mean 
field dynamo problem has focused on unit Prandtl number Pr = A/77 (the 
ratio of the magnetic to kinetic viscosities.) In nature, this number is rarely 
near unity. In the Sun Pr « 1, whereas in the Galaxy Pr >> 1. There are 
currently several numerical [19,42,43] and analytic studies [43] addressing the 
role of the magnetic Prandtl number. For unit Prandtl number, the magnetic 
and kinetic energies approach equipartition from the the input scale down to 
the dissipation scale. For Pr » 1, Ref [41] finds that the magnetic energy 
dominates the turbulent energy at scales below the viscous cutoff scale, and 
that the magnetic energy initially peaks on these tiny scales. With time, 
the peak seems to migrate back toward the input scale, although runs have 
not been done long enough to see what the saturated state is like. These 
simulations do not have a large enough inertial range, and so are far from 
being able to assess the power build up at scales even larger than the forcing 
scales, as required to test the MFD. nor do they consider helical forcing. 

Ref [21] finds that the non-local inverse cascade is not much affected by 
the large Pr when compared to the small Pr, but intially sees a similar 
development of the small scale magnetic energy as in Ref [42]. The dynamic 
range is also limited. Refs [43] appears consistent with this. The work also 
suggests that the contribution to a from the current helicity might exceed 
that from the kinetic helicity for large Prandtl number. More work is needed 
to understand the effect of large and small Prandtl numbers. 

V. CONCLUSIONS 

We have emphasized that the MFD represents a framework for under- 
standing an inverse cascade of magnetic helicity: Kinetic helicity imposed 
at small scales pumps magnetic helicity from small to large scales through 
an inverse cascade [21,22,44]. This process seems to be non-local [22] in 
that helicity "jumps" directly from small to large scales. Accompanying the 
magnetic helicity is a growth of magnetic energy. 

We suggest that the MFD can in principle proceed much more rapidly for 
real astrophysical rotators, when compared to simulations in periodic boxes 
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in which the turbulence is homogeneous. For such simulations, surface terms 
are ignored and so the magnetic helicity is nearly conserved for large magnetic 
Reynolds numbers. Generating magnetic helicity and magnetic energy on the 
large scales from an MFD then requires generating a compensating magnetic 
helicity at or below the input scale. In a steady state, the small scale helicity 
would drain only through resistive dissipation. In this case the growth rate 
of the large scale field is resistively limited. 

But in real astrophysical systems, the boundary terms relax the helicity 
conservation constraint. The growth rate of the large scale helicity (and 
thus large scale field) is then limited by the rate at which the compensating 
helicity flows out the boundary. This rate can in principle be much faster 
than the resistively limited rate. The message is that the boundary terms are 
likely important for astrophysical MFD's and so backreaction studies which 
show suppression but do not allow boundary terms leave the ambiguity as 
to whether the suppression is actually due to the dynamic backreaction or is 
simply due to the boundary conditions. 

Note that the role of shear for the generation of magnetic 

A steady flow of magnetic helicity into the corona is expected for an astro- 
physical rotator harboring a vigorous a — Q mean field dynamo. The helicity 
escape rate leads to a lower limit on the total magnetic energy deposition 
into the corona. When the corona itself is turbulent, there should also be 
an inverse cascade of magnetic helicity in any wind driven outward, thus the 
dominant magnetic helicity scale would appear to increase on increasingly 
large distances from the source. 

In addition to the conclusions above, we discussed a number of open 
issues: the choice of boundary conditions for simulations, the need to under- 
stand the boundary diffusion in real systems, some implications of previous 
analytic work, the predicted variability for a corona, and the role of the 
magnetic Prandtl number. 
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